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Abstract We provide a brief overview of the modelling of the atmospheres of red 
giant stars with the 3D radiation-hydrodynamics code C05B0LD. We emphasize 
aspects where 3D modelling provides additional insight beyond standard hydrostatic 
ID models, and comment on present modelling challenges. 



1 Introduction 



Modelling of three-d imensional (3D ) atmo spheres of cool stars is an active field 
of development (e.g. Nagendra et al. . 20091) . and particularly 3D models of atmo- 
spheres of red giant (RG) stars are just on the verge of bec oming available for ap- 
plication to astrophysical problems. In an early application. iKucinskas et al J (120051) 
used a 3D RG mod el to estimate color corrections due to thermal inhomogeneities; 



Collet et al.l (120071) considered a set of eight giant models to investig a te the impact 
on line formation and abundances. More recentlv. iFrevtag & Hofner (2008) devel- 
oped model atmospheres of AGB stars and their winds. iDupret et al..(i2009i) derived 



the en ergy input to solar-like oscillations in giants from 3D models, iRamfrez et al 



(1201 Ol) studied convective li ne-shifts in the metal-p oor RG HD 122563 and com 
pared them to a 3D model, IChiavassa et al.l (1201 lb applied global 3D models to 



assess effects of photometric and related astrometric variability, and iPasquini et al 



(1201 11) took recourse to 3D dwarf and RG models to correct for convective blueshifts 
in high-precision, spectroscopic radial velocity measurements. While fairly exhaus- 
tive, the list of examples is still quite short, but illustrates already the variety of 
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possible applications of 3D RG models. At the moment efforts are under way to 
cover the Hertzsprung- Russell diagram with 3D model atmospheres inclu ding stars 



in the red-giant branch (ILudwig et al.Ll2009atlTrampedach & Steinll201 lb 



In the following we are going to focus on aspects related to abundances from 
3D models, and the theoretical calibration of the mixing-length parameter ttMLX 
from 3D model atmospheres. Here, we are not so much presenting results as rather 
pointing out problems which are still lingering. We finally add some comments 
about predictions of the photometric micro-variability which are of interest in the 
context of high-precision photometry missions like CoRoT. 



2 Our model atmosphere codes 

2.1 3D model atmospheres - C05B0LD 



Our 3D mod el atmospheres were calculated with the radiation-hydrodynamics code 
C05B0LD dFrevtag et al.ll2002l:IWedemever et al.ll2004l:lFrevtag et al.ll201 lb . The 
code solves the time-dependent equations of compressible hydrodynamics coupled 
to radiative transfer in a constant gravity field in a Cartesian computational domain 
which is representative of a volume located at the stellar surface. The equation of 
state takes into consideration the ionization of hydrogen and helium, as well as the 
formation of H? molecules according to Saha-Boltzmann statistics. Relevant ther- 
modynamic quantities - in particular gas pressure and temperature - are tabulated 
as a function of gas density and internal energy. The multi-group opacities used 
by C05B0LD are based on monochromatic opacitie s stemming from the MARCS 



stellar atmosphere package dGustafsson et al.L 120081) provided as function of gas 



pressure and temperature with high wavelength resolution. The opacities have been 
calcul ated assuming solar elemental abundances according to iGrevesse & Sauval 
(Il998h. with the e xception of CNO for which values close to the recommendation 
of lAsplundl ( l2005b are adopted (specifically, A(C)=8.41, A(N)=7.80, A(0)=8.67). 
The metal abundances were scaled according to overall metallicity of the model as- 
suming an enhancement of the a-elements by H-0.4 dex at metallicities [M/H] < — 1 . 

In our RG models we typically use a number of 140 x 140 x 150 to 160 x 160 x 
200 points for the hydrodynamical grid. The decision about the resolution primarily 
hinges on the effective temperature of the model, hotter models usually require a 
higher resolution. The wavelength dependence of the radiation field is represented 
by 5 multi-group bins in the case of solar metal licity, and 6 bins at sub-sola r metal - 
licities, following the procedures l aid ou t by Nordlund (1982); Ludwid (1992); 
Ludwig et al.l ( 119941) ; IVogler et al.l ( 120041) . For test purposes we have calculated a 
few models with more bins. Since it is of relevance for the discussion later, we 
emphasize that all opacity sources - including scattering opacities - are treated as 
true absorption. The sorting into wavelength groups is done applying thresholds in 
logarithmic Rosseland optical depth {+oo,0.0, — 1.5, —3.0, —4.5, — oo} forthe 5-bin, 
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and {+00,0.1,0.0, —1.0, —2.0, —3.0, — oo} for the 6-bin schemes. In all but one bin a 
switching between Rosseland and Planck averages is performed at a band-averaged 
Rosseland optical depth of 0.35; in the bin gathering the largest line opacities, the 
Rosseland mean opacity is used throughout. The decisions about number of bins, 
and sorting thresholds are motivated by comparing radiative fluxes and heating rates 
obtained by the binned opacities in comparison to the case of high wavelength res- 
olution. 



2.2 ID stellar atmospheres - LHD 



Due to still present Umitations in the realism (e.g. by the limited wavelength resolu- 
tion) of 3D model atmospheres it is often advantageous to work differentially, and 
express 3D effects relative to a ID comparison structure. To this end we developed 
a ID stellar atmosphere code called LHD which employs the same opacities and 
equation-of-state as the 3D code C05B0LD. The convective en ergy transport i s 
modelled in the framework of mixing-length theory as described in lMihalasI (119781) . 
The result ing ID stratification s are in hydrostatic and radiative-convective equilib- 
rium. See ICaffau et al.l (120071) for more details on our approach of deriving abun- 
dance corrections. 



3 RG abundances and the issue of scattering 



Collet et al.l (l2007h presented 3D- ID abundance corrections for RG models at effec- 
tive temperatures of around 5000 K, \ogg = 2.2, and meta llicities ranging from sola r 
to [M/H] = -3, using the 3D code o f Nordlund & Stein (IStein & Nordlund[ll998l). 
Two si milar studies were presented by Dobrovolskas et al. ( 20 id) and Ivanauskas et al. 
(I2OIOI) who used C05B0LD and LHD models at about T^g = 5000 K, logg = 
2.5, with metallicities down to [M/H] = —3 to derive 3D-1D abundance correc- 
tions. While the two later studies confirm the results of Collet and collaborators, 
showing that generally the magnitude of 3D- ID abundance corrections becomes 
larger towards lower metallicity, the quantitative agreement is not satisfactory. The 
C05B0LD-based abundance corrections are usually noticeably smaller in magni- 
tude, in particular at the lowest metallicities. Obviously, this is an unfortunate situ- 
ation, and one would like to see a higher degree of consistency among results from 
different 3D codes. 



In a recent paper. ICoUet et al.l (1201 Ih suggested the treatment of scattering in the 
simulations as the reason for the discrepant abundance corrections for RGs at low 
metallicity. The main scattering process is Rayleigh scattering by neutral hydrogen. 
This is perhaps the simplest case of scattering and can be modelled as coherent 
isotropic scattering in the continuum. Collet and collaborators implemented a proper 
treatment of this kind of scattering in 3D. They also put forward an approximate 
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treatment of scattering by simply leaving out scattering contributions in the binned 
opacities in the optically thin regions. They showed that this approximate treatment 
provides results in close agreement with the exact treatment. They further performed 
a comparison with the case where scattering is treated as true absorption - as is 
the case in the C05B0LD models. Their models of 2007 used the approximate 
treatment of scattering. The models show a sensitive dependence of the resulting 
temperature stratification on the treatment of scattering. In their [M/H] = —3 RG 
model, the difference amounts to 600 K at optical depth log Tross = —4 in the sense 
that a proper treatment of scattering leads to cooler structures in comparison to 
treating scattering opacities as true absorption. 
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Fig. 1 Comparison of ttie mean temperature structures of two different 3D C05B0LD hydrody- 
namical model atmospheres (solid) and associated ID LHD models (dashed), computed with a 
different treatment of radiative transfer. In the first case (dark [red] curves), the continuum scatter- 
ing opacity is treated as true absorption opacity, while in the second case (light [green] curves), the 
continuum scattering opacity is ignored in the optically thin layers. For the 3D models, averaging 
was performed on surfaces of constant Rosseland optical depth and over 70 equidistant snapshots 
covering a total of 140000 s. 



It appears plausible that the differing 3D- ID abundance corrections are a conse- 
quence of the different thermal structures resulting from the different treatment of 
scattering in the C05B0LD and Nordlund-Stein class of models. To test this idea, 
we calculated a RG model with the same atmospheric par ameters as before bu t with 
the approximate treatment of scattering as suggested by ICollet et alj (1201 In . Fig- 
ures [T| and |2] illustrate the outcome. The most striking aspect is that our models 
show a very much reduced sensitivity to the treatment of scattering in comparison 
to the models of Collet and co-workers. The approximate treatment of scattering 



Page: 4 job: ludwig_roma macro : svmult . els date/time : 24-Aug-2011/0 : 16 



3D Model Atmospheres of Red Giant Stars 



1000 



800 



600 



400 



200 



d3t50g25mm30n01, d3t50g25mm30n02 



true absorption 
approx. scattering 



-5 -4 -3 -2 -1 

log Xross 

Fig. 2 Total rms temperature fluctuation AT^msxat as a function of Rosseland optical depth for 



the two 3D models shown in Fig.[T] computed as ^Trms.tot = y {T^)x,yj — {T)-^y,, where 
denotes horizontal averaging over surfaces of constant Rosseland optical depth and over time. We 
have verified that the amplitude of the total temperature fluctuation is completely dominated by the 



spatial temperature variations: AT^msMt ~ 



{T)ly 



leads to a structure which is only 120 K cooler at log Tross = —4, in comparison 
to w 600K found by ICollet et alj (1201 ih . This also carries over to the temperature 
fluctuations which are little affected by the treatment of scattering (see Fig.[2]i. We 
already emphasized the importance of a differential approach, and Fig.[T]also shows 
the effects on the associated ID LHD models. Temperature differences between the 
ID and 3D models at given optical depth are changing even less. While we did 
not perform spectrum synthesis calculations yet to derive new abundance correc- 
tions, we consider it as unlikely that the modest changes in the thermal structure 
can change our abu ndance corrections so much that they become consistent with 
the values of lCoUet et al. (2007i) . 

The situation remains puzzling. The very different sensitivity to the treatment of 
scattering is difficult to explain. We only can hint at the differences in the calcula- 
tion of the band-averaged opacities in the various codes: Collet and collaborators 
use intensity-averaged opacities in the optically thin regions, while we use Planck- 
averages - except for the band collecting the strongest lines where a Rosseland 
average is used throughout. We speculate that these choices, together with the defi- 
nition of the opacity bins, may have a significant influence on the resulting thermal 
structures and their sensitivity to the treatment of scattering. 
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4 The calibration of cumlt and turbulent pressure 

It is well known from the theory of stellar structure that convection is generally 
an efficient means of transporting energy, and that it establishes a thermal struc- 
ture close to adiabatic. Only in the vicinity of the boundaries of convective regions 
noticeable deviations from adiabaticity occur In convective envelopes of late-type 
stars the upper boundary of the convective envelope - usually located close to or 
even in the optically thin layers - constitutes the bottle-neck for the energy trans- 
port through the stellar envelope assigning a special role to it. Despite its small 
geometrical extent and low mass, it largely determines the properties of the convec- 
tive envelope as a whole. It is the value of the entropy of the adiabatically stratified 
bulk of the convective Senv which is most important from the point of view of stellar 
structure since it influences the resulting radius and effective temperature of a stellar 
model. Senv is controlled by the efficiency of convective and radiative energy trans- 
port in the thin, superadiabatically stratified surface layers. 3D model atmospheres 
can be applied to model this region, and allow to quantify the mutual efficiency of 
the convective and radiative energy transport, and to predict Senv Comparing the 
model predictions to standard ID models based on mixing-length theory (MLT) the 
value of Ssnv can be trans l ated into a corresponding mixing-length parameter aMLT 
dTrampedach et al.lll999l:lLudwig et al.lll999ll2008l) . 



In stellar evolution calculations the free mixing-length parameter is usually cal- 
ibrated against the Sun. However, it is unclear whether mixing-length theory pro- 
vides a suitable scaling of the convective efficiency at constant aMcr across the 
Hertzsprung-Russell diagram. The depth of the surface convective envelope and 
the related OCmlt can be constraint by asteroseismology. However, degeneracies 
with other parame ters often make it difficult to obtain a unique solution (e.g. 
Goupil et al.L [201 lb . Hence, it would be useful to have an independent estimate 
available which 3D models can provide in principle. 

In main-sequence models turbulent pressure plays generally only a minor role 
but becomes relatively more important towards lower gravities - and causes trou- 
ble when one is interested in a well-defined calibration of the mixing-length pa- 
rameter. Figure [5] shows the average temperature profile of a 3D red giant model 
(Teff ~ 3600 K, logg=1.0, [M/H]=0.0) in comparison to standard ID model atmo- 
spheres of the same atmospheric parameters. While the turbulent pressure Pturb is 
naturally included in the 3D models, it is modelled in a ad-hoc fashion in ID mod- 
els, assuming a parameterisation Pturb = ftmbPv^, where /tuit is a free parameter of 
order unity, p the mass density and Vq the convective velocity according MLT. 

Figure [3] shows that it is essentially impossible to reproduce the mean thermal 
profile of the 3D model with any of the ID models - irrespective of the choice 
of aMLT and /turb- The failure is related to the local nature of MLT confining the 
action of the turbulent pressure gradients strictly to the convectively unstable re- 
gions. While formally one can still match Senv of the 3D model by a ID profile 
with suitably chosen aMLT and/or /turb such a match becomes physically little mo- 
tivated, and is unlikely to provide a robust scaling with changing atmospheric pa- 
rameters. An improved ID convection description including non-local effects like 
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Fig. 3 Entropy-optical depth profiles (horizontal and temporal average) of a 3D red giant model 
(thick solid line) in comparison to ID stellar atmosphere models of different Kmlt leaving out (thin 
solid lines) or including (dashed and dashed-dotted lines) turbulent pressure. The lines are labeled 
with values «MLT//tuib (details see text). The horizontal dashed line indicates the value of Senv as 
predicted by the 3D model. 



overshooting is clearly desirable to handle this situation. Empirical calibrations of 
OCmlt using giants are likely to suffer from ambiguities related to the way turbulent 
pressure is treated in the ID models. One may take the result as an indication that 
taking recourse to ID models is not warranted, and one may give up the benefits 
of a differential approach by relating 3D to ID structures. Alternatively, one may 
take the absolute entropy of the convective envelope (perhaps translated to equiv- 
alent pressure-temperature pairs) as predicted by the 3D model as constraint to be 
matched in ID stellar structure models. 



5 Granulation-related photometric micro- variability 



High-precision photometry of satellite missions (foremost Most, CoRoT, Kepler) 
allow the detection of stellar variability associated with the random changes of the 
granulation pattern on the surfaces of late-tape stars - by asteroseismologists usu- 
ally referred to as "granular background noise". 3D model atmospheres represent the 
granulatio n pattern in detail and allow to predict the power spectrum of the var iabil- 



ity signal (ITrampedach et al.Lll998l : ISvensson & Ludwia , l2005l:lLudwia,l2006l) . De 



spite this possibility, no comprehensive theoretical study has been conducted so far 
One of the reasons is that long time series need to be calculated to collect sufficient 
statistics, which is computationally demanding. The F-dwarf HD 49933 - a promi- 
nent CoRoT-target - is an exception for which iLudwig et al. I (l2009bh performed 
an analysis. However, the growing body of observational data in particular for giant 
stars s hould motivate further efforts in this direction. Recentlv. lKieldsen & Bedding 
(1201 In suggested a new scaling relation for the amplitude of solar-like oscillations, 
and also discuss the scaling of the granulation background signal. It would be in- 
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teresting to see whether 3D model atmospheres can lend further support to the sug- 
gested relations. 

To illustrate the feasibility, we show in Fig.|4]a rough comparison of the photo- 
metric variability between the RG HD 181907 (HR7349) and predictions from two 
3D models. The plot focuses on the frequency region where the granulation-related 
signal is expected. CoRoT acquired a high-quality time series for HD 181907; 
Carrier et all tOldfi give atmospheric parameters 4780 ± 80/2.78 ± 0. 16/-0.08 ± 
0. 10 (7;ff/logg/[M/H]). The two 3D models have atmospheric parameters 4500/2.5/0.0 
and 5000/2.5/0.0, bracketing the star in effective temperature, as well as having 
comparable surface gravity and metallicity. Although no dedicated modelling was 
performed, the spectra appear quite similar. 




0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 

Frequency v [mHz] 

Fig. 4 Comparison of power spectra of photometric vaiiability. Black line: HD 181907 as observed 
by CoRoT. The predicted power spectra from two 3D models are depicted with light [green ] 
(Jeif = 4500 K) and dai'k grey [red] (Feff = 5000 K) lines. Further details see text. 



6 Concluding remarks 

3D model atmospheres of cool stars, including red giants, have reached a level of 
realism which allows a direct confrontation with observations. In some areas they 
allow to make predictions beyond the capabilities of classical ID models. However, 
as we have seen there still exist modelling challenges, and last but not least quite 
some work is still necessary to fully exploit the potential of such models. 
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